Knot solitons in modified Ginzburg-Landau model 
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We study a modified version of the Ginzburg-Landau model suggested by Ward and show that 
Hopfions exist in it as stable static solutions, for values of the Hopf invariant up to at least 7. We also 
find that their properties closely follow those of their counterparts in the Faddeev-Skyrme model. 
Finally, we lend support to Babaev's conjecture that longer core lengths yield more stable solitons 
and propose a possible mechanism for constructing Hopfions in pure Ginzburg-Landau model. 
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I. INTRODUCTION 

Topological solitons have long enjoyed widespread in- 
terest within many fields of physics, including such seem- 
ingly distant subjects as cosmology, condensed matter 
physics and particle physics. Of these, the most com- 
mon example are perhaps the Abrikosov vortices in type 
II superconductors. Therefore, it is crucial to under- 
stand the properties and existence of topological solitons. 
The purpose of this work is to provide further informa- 
tion about their presence in Ward's modified Ginzburg- 
Landau model. 

While many topological solitons are pointlikc or two- 
dimensional, there are also extended three-dimensional 
topological solitons; one particular class of them is called 
Hopfions, their name arising from the topological invari- 
ant associated with them, called the Hopf invariant. The 
archetype model supporting topologically stable closed 
vortices, is the Faddeev-Skyrme model [ll-Q. While this 
has been know for quite some time, it has recently be- 
come more relevant with the discovery of topological 
insulators and the possibility of the existence of Hop- 
fions in them It is also known [ll|, [T^] that that 
the Faddeev-Skyrme model can be embedded into the 
Ginzburg-Landau model, giving rise to a conjecture in 
[l3 | that the two-component Ginzburg-Landau model 
could, due to this embedding, also support the same 
topological structures as the FS model does. 

Previous research has not been able to reach a conclu- 
sion on the stability of Hopfions in the Ginzburg-Landau 
model. There has been one positive result in a restricted 
model [l3| but it has not been confirmed in other inves- 
tigations using the full model (T^ - [l6| . This suggests that 
any Hopfions in the Ginzburg-Landau model are very 
hard to find, so studying closely related models is not 
only relevant in itself, but may also provide clues as to 
how to find Hopfion in the Ginzburg-Landau model. In 



the modified Ginzburg-Landau model, however, stable 
Hopfions have been discovered [iB] , but these works 
have only explored the stability of Hopfions at Hopf in- 
variant one. We will show that Hopfions exist as local 
minimal energy configurations in the modified model up 
to at least the value 7 of the Hopf invariant and that 
these share several important features with their coun- 
terparts in the Faddeev-Skyrme model. Finally, we pro- 
vide support for a conjecture by Babaev that solitons 
with longer cores are more stable and propose a possible 
way to construct Hopfions in the unmodified Ginzburg- 
Landau model. 



II. THE MODEL 

The static Abelian Higgs model with two charged 
Higgs bosons has the same mathematical form as the 
Ginzburg-Landau model with two flavors of Cooper pairs 
or superfiuids. We will use the following notations. The 
indices run as follows: j,k,l e {1,2,3}, a £ {1,2}, 
and the fields are ^I* = {ipi ■02)^, Fjk = djAk — dkAj, 
~ e^^^dkAi and the gauge-covariant derivative has 
the form Dj = dj — igAj; for short, we will also write 
D = V — igA. With these notations, the standard static 
energy density of the two-component Ginzburg-Landau 
model can be written as 



(1) 



where we have used SI units. The exact form of the 
potential depends on the physical context, but for the 
purposes of this article, all that is required is that it 
maintains the SU (2) symmetry of the model and enforces 
the condition ||4'|| = constant ^ at some limit of the 
parameters (now 77) of the potential. Here we have used 
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The embedding of Babaev et al. [T^l is useful to 
demonstrate how closed (or knotted) vortices might ex- 
ist in Ginzburg-Landau model. These will be defined in 
terms of the fields tpa, leaving the gauge field A free. 
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The embedding requires that ||^'|| > everywhere, but 
his is not enough to reveal the closed vortices and further 
conditions need to be met as we shall next describe. In 
all the topological considerations that follow we assume 
that 1 1^*1 1 > everywhere and 'i/ can be thought of as 
normalized by ^ — > \E'/||4'||, which implies that v]/ e S^. 

Maps S'^ — >■ CP^ = fall into disjoint homotopy 
classes, the elements of 7r3(S^) and it is well known that 
7r3(S^) = Z. Thus, for each map — > CP^ we can 
assign an integer, called the Hopf invariant, which tells 
us which element of 7r3(S^) that belongs to. Similarly, 
maps S'^ — ?> S'^ belong to elements of 7r3(S^) = Z; this 
integer is called the degree. Suppose one has two maps: 
* : ^ and the Hopf map h : ^ CP^ ^ S^. Then 
it can be shown that the Hopf invariant of ft. o v]/ equals 
the degree of Vl/. Since we work with maps 5^ — > S"^, the 
relevant topological invariant is the degree, but there is 
always an associated Hopf invariant, iJ, which has the 
same value. Topological solitons, i.e., stable static solu- 
tions with an associated Hopf invariant, are called 
Hopfions and also knot solitons, due to their general 
shape at higher values of H. Next, we will see how it is 
the presence of the Hopf invariant, not the degree, which 
gives rise to closed vortices and knot solitons. 

Without loss of generality, we can choose = 
(1, 0). Now, the Hopf map takes ^ 4>^ =(0, 0, 1) € 
and we define the soliton core as the preimage 

It is natural to ask whether one of these Hopfions is the 
global energy minimum for each H as is the case in the 
Faddeev-Skyrme model. The answer to the question is, in 
general, negative. The fact that A is left free, means that 
there is no nontrivial topology imposed on it and since 
in the vacuum of the Ginzburg-Landau model A is pure 
gauge, the magnetic field energy can vanish in all cases. 
This in turn means that there is no longer a fourth-order 
derivative in the energy density ([1]) and Derrick's the- 
orem [l^ states that no stable, topologically nontrivial 
solutions of the field equations with nonzero energy exist. 
This fact has been observed many times, by various au- 
thors, including [13, [SHO]- It has also been seen in nu- 
merical work, that the magnetic field energy indeed does 
vanish and with it, the soliton itself [Tjilal- In short, 
the knot soliton can always be undone by decreasing the 
magnetic field to zero and radially shrinking it. However, 
if the starting point is chosen suitably, this process may 
involve temporarily increasing the energy of the configu- 
ration. Thus, local minima may still exist. Indeed, they 
were shown to exist in the modified model in [IJ, [l^ , 
while in [l^ the collapse of the magnetic field described 
above was observed and no stable solutions were found 
in pure Ginzburg-Landau model. 

The crucial ingredient to finding these local minima 
is to prevent the collapse of the magnetic field in order 
to obtain stable topologically nontrivial configurations 
in the model. There are physical arguments that suggest 
there might exist physical processes that prevent the col- 
lapse, but here we follow the path set out by [ij], where 



the energy density is modified by adding the term £w 
(we denote '^'^ = (-01 tl'2))' 



£glw ^ \\\Dn^ + \\\^ -K A\\ 



=£p 
^(^1,^2) 



(3) 



Denoting for any subscript z: = Jd^xSz we finally 
have the total energy 
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The extra term, when the parameters K,r] ^ 00, ensures 
that the model becomes exactly the Faddeev-Skyrme 
model 



(5) 



where is a normalized real valued field ^ S^. There- 
fore the model supports Hopfions, at least asymptotically. 
This limit of k — > 00 was apparently first observed by 
Hindmarsh albeit in a slightly different context. 



III. NUMERICAL METHODS 

We will now turn to our numerical investigation. All 
the computations are done using the methods and pro- 
grams described in [l^. In brief, we have used gauge- 
invariant discretization of the energy density with sim- 
ple forward differencing scheme for derivatives. We then 
compute the gradient of the discrete energy density with 
respect to the fields and use a conjugate gradient 

algorithm to find a local minimum from a given ini- 
tial configuration. The initial configurations of ^ for 
lower values of H are constructed using the methods of 
[l5| . that is, denoting the coordinates of by Xi and 
= a;? + x?, + x^, we set 
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+4xi 



(6a) 
(6b) 



These satisfy > everywhere and indeed = 1 
everywhere. The degree of this configuration is deg ^' = 
pq, so we can easily construct configurations of any given 
degree - or Hopf invariant - when one thinks about the 
embedding. For H > 5, the initial configurations de- 
scribed in section 3 of [2lj were used. 

We also need an initial configuration for A, which we 
take to be always 

Ai{x) = -lni^-V'i(x) (i'i{x + hii) - tpi{x)^ 

-^p2{x)(jj2{x + hei) -ip2{x)U/h, (7) 
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where h is the lattice constant. The choice of initial value 
for A is largely irrelevant since it is has the same homo- 
topy (trivial) in any case, but Eq. dH) convergences much 
faster relative to the obvious initial configuration Ai = 0. 

In order to save computer resources (time), we reused 
previously found minimum configurations to look for new 
ones: once a stable minimum was found for some values 
of iJ, K, r], we sometimes used that solution as an initial 
configuration for some new values oi K,ri (of course, we 
cannot change H in this manner). This initial configu- 
ration converges faster than a fresh configuration set up 
using Eq. (O , which we interpret to indicate that the old 
minimizer is in some sense "closer" to the new solution 
than the fresh configuration. Furthermore, it is often 
true, that such a reused minimizer converges to a stable 
solution even when the fresh initial configuration does 
not, that is, the fresh configuration is sometimes outside 
the attraction basin altogether. This emphasizes the fact 
that we can never be certain that there is no stable solu- 
tion for a given pair k, r/ using our methods: the attrac- 
tion basin could simply be so tiny that constructing an 
initial configuration within it is nearly impossible with- 
out some additional knowledge about its location within 
the configuration space. We will, however, present strong 
evidence pointing to the nonexistence of stable solutions 
for certain values oi H, n^rj. 

In modern numerical work, data post-processing is in 
an increasingly complex role. We will now describe the 
data post-processing used in this work. All the post- 
processing is done using the MayaVi Data Visualizer 
[22| . with some extensions we have implemented our- 
selves in the python language. Even though MayaVi 
is designed for visualization, it contains routines gen- 
erally useful for post-processing, such as interpolation, 
which we use heavily. For an isosurface plot we compute 
(j) = ho^i and then simply ask MayaVi to find, using in- 
terpolation, the surface satisfying the conditions for the 
isosurface in question. 

Computing the core length is a more complex opera- 
tion since interpolation cannot be directly used to find 
the minimum values of data. However, we can work 
around this limitation using the fact that \\(f)\\ = 1. We 
use MayaVi's routines to first filter out points where 
(/)3 > 0, since we are not interested in these. Then we 
apply MayaVi's contour finding routine two successive 
times to the remaining data: first, we find the contour 
4>i = 0, i.e. points where 02 + = 1 ^"^^ then from 
this data, we find the contour 02 = 0, i.e. points where 
0| = 1, but since we excluded points where 03 > 0, what 
is left is the core 03 = —1. MayaVi represents this as a 
polygon, whose circumference we then compute to give 
the core length. The main advantage of this method is 
that MayaVi is able to interpolate the contours at each 
stage, giving smoother cores than would be possible with 
direct methods. The accuracy of this procedure is very 
good: at a lattice with 60^ points and lattice constant 
h — 0.1, the interpolated core length has an error of less 
than 0.03%. 



As an example of previously described routine, we 
present the trefoil knot shaped local minimum at iJ = 7, 
K = /y = 10 in Fig. [TJ The other preimage is constructed 
in a similar fashion. 




Figure 1: The isosurface of 03 = —0.5 for the local 
energy minimum at k ~ rj — 10. The solid curve inside 

the isosurface is the soliton core and the other solid 
curve is another preimage, namely 0~^(— 0.87, 0, —0.5). 



IV. RESULTS 

We have investigated the properties of stable solutions 
in the modified Ginzburg-Landau model for H G [1,7]. 
All numerical results have been obtained in a standard- 
ized lattice with 360'^ lattice points and lattice constant 
h = 1/12. Since this is different from the lattice of (l6j . 
we also recomputed the H = 1 case in the new standard 
lattice in order to be able to compare the results for all 
values of H. Some checks were made in lattices of 480"^, 
h = 1/12 and 180^ /i = 1/6 as wcU. 

We chose to investigate the solutions (stable or other- 
wise) for each H along two lines: k = 10 and 77 = 10. 
As expected, for both lines, we find that there is a lim- 
iting value of K (rj) below which no stable solutions can 
be found. These values depend on H. For H = 1 our 
results arc in agreement with [3, HB] and stable solutions 
are also found for all investigated H > 1. 

This procedure gives us two points on the sta- 
ble/unstable boundary investigated in [l6j . and also in- 
formation about how the stable solutions change with 
changing k and 77, which is not so easy to discern from 
solutions following the actual stable/unstable boundary, 
where both k and 77 are changing. 



A. Error considerations 

Derrick's theorem [31 implies that for any static solu- 
tion of the field equations which is stable against uniform 
scaling of space, a virial theorem must hold. In this case, 
it takes the form Eb = Ek + Ew + SEp. Obviously, 
for a numerical approximation, the equation will only be 
satisfied to within some tolerance. This tolerance would 
ideally be deduced from the discretization errors, effects 
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of a finite computational domain and the tolerance used 
to determine convergence, but this seems an insurmount- 
able task. Instead, we have computed the same solution 
with different lattices (varying both size and density) to 
give us an estimate of the accuracy. This then gives us 
an idea of the tolerance in the virial theorem which is 
achievable at a given lattice. We use this tolerance to 
estimate the errors in the results. 

Apart from error estimates for each simulation, it is 
important to note what this procedure tells about our 
results in general. Fig. [2] depicts the energies and core 
lengths of the H = 7 simulations for k = 10 (panel [^a|) 
and for r] = 10 (panel [2b|) with two difi'erent lattice sizes. 
As can be seen, the differences are minimal apart from 
the lowest values of rj, where the energy drops signifi- 
cantly faster in the smaller lattice. The difference in the 
soliton itself is not easily seen until r/ = 0.05, where we 
only find a stable soliton in the larger lattice. This is due 
to the boundaries exerting pressure on the growing (as rj 
decreases) soliton, thus destabilizing it. 




IQ ' 10" 10' 10^ 10- 

(a) K = 10 



(b) r? = 10 



Figure 2: Total energy and its error limits in normal 
size and large grid for charge Q = 7 trefoil. Triangles 
pointing down (solid blue curve) represents 360'^ and 
triangles pointing up (solid green) 480'^ lattice, dashed 
(blue) line the core length in first case and dash-dotted 
(green) in second case. Energy values are on the left 
and core length values on the right y-axis. 

The same effect is present at increasing k, but here it 
cannot destabilize the soliton since increasing k takes us 
towards a more Faddeev-Skyrme -like system, where the 
soliton is stable. The effect, however, is large enough to 
give rise to fast decrease in accuracy. This decrease is 
the reason for the disparity on the ranges covered: rj G 
[0.05, 1000] and k e [6, 14] in this work. 

As noted in Section IIIH the errors in the core length 
estimates are negligible as long as the core is relatively 
smooth. 



B. Soliton energy 

Let us now describe in detail the effects of varying the 
two parameters on the soliton solutions. First, it should 



be noted that the energies of the solitons follow the same 
E oc H^/^ as those of the Faddeev-Skyrme model. There- 
fore, in what follows, we always normalize the energy: 
E -> E/H^/\ 

Because of the asymptotic limit ([5]), one expects the 
energies of the solitons to approach the limit set by the 
solitons of the Faddeev-Skyrme model as k, ?/ — > oo. This 
expectation is found to be true: for each value of the Hopf 
invariant at small values of k or ry, the energy is well be- 
low the limit and starts growing towards the value of the 
Faddeev-Skyrme model when either of the parameters is 
increased, as seen in Figure [31 where we plot the ener- 
gies of all the identified local minima for all investigated 
values of as a function of (panel [5a)) and k (panel 




10-' 10" 10' lO-' 10- 



(a) K = 10 



(b) = 10 



Figure 3: Plots of the normalized energies [E/H^/'^) of 
the minimizers for each H = 1 (blue circle), H = 2 
(green triangle) , H = 3 (red square) , H — A (cyan 
diamond) and H ^ 7 (magenta pentagon). We also 

include a single data point for H = 5 (yellow cross) and 
H = 6 (black plus). 

It is worth noting at this point that our program seems 
to systematically slightly overestimate the energy: for 
H = 4, our values at k = 14,7^ = 10 and k = 10, ?y = 
1000 are about 5% above those reported for the Faddeev- 
Skyrme model where an underestimate of about 1% 
was reported. We believe that since our results show 
a clear approach towards an asymptotic value, this is 
an acceptable accuracy for such a simple differentiation 
scheme as used here. The accuracy could be improved 
by simply using smaller lattice spacings, but that seems 
unnecessary for the purposes of this work and would be 
computationally too expensive for such a large number 
of simulations. 

We also note that at values ofry = 10, k > 10 and 
T] > 10, K = 10, the order of the normalized energies 
is the same as in the Faddeev-Skyrme model [11-0, HJ , 
further emphasizing the close relationship of the models. 
It is interesting to note, that the same order also appears 
in the extended Faddeev-Skyrme model [23[ . We include 
two extra data points at k = 77 = 10 for H = 5,6 to 
further demonstrate this. No other simulations were done 
at iJ = 5, 6. At lower values of the rj, the order changes, 
but simultaneously the error bars grow significantly. For 
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small K the order would seem to persist, but this is not 
the case: the boundary of the region of stable solitons is 
reached at different values of k for different values of H; 
for example, as H = 1 reaches the boundary before H = 
2 it will necessarily have lower energy below some limiting 
value of K, thus disrupting the order as is explicitly seen 
to occur for small ?/. 

We want to emphasize the fact that the boundary we 
discover is a bound on the values of k, rj above which sta- 
ble solitons can be found. There is no evidence of their 
existence below the boundary, but it cannot be ruled 
out using our methods. Also, even for our methods, 
the bound can be pushed slightly downwards by using 
more accurate lattices, but instability still eventually oc- 
curs (up to what is computationally feasible) as seen for 
if = 1 in 0. 



C. Soliton core length 

Next, we turn to the length of the soliton core, i.e. the 
length of the curve (/)-i(0, 0, -1) S S3=K3u{oo}. 

Sutcliffe found that core length for Faddeev-Skyrme 
Hopfions follow curve ^H^^'^ where his data was fitted 
to the curve to produce value 7 = 7.86 [2l[. We plot 
the core lengths of our soliton solutions in Figure |H It 
is immediately obvious from panel I4bl that the solitons 
collapse to zero size as k — >■ and as k grows, the soliton 
sizes approach some asymptotic values.. The behavior of 
the core length is more complex in panel I4al where k is 
held constant. Here the core lengths increase, seemingly 
without limit, as ij decreases, but do not seem to be- 
have monotonically. However, even though the method 
used to determine the core length from the computational 
data is extremely accurate, this gives no information as 
to how accurately the numerical core approximates its 
continuum counterpart. The seemingly increasing size 
of the core for large values of rj falls within the esti- 
mated accuracy of the program, so it can very well be 
a numerical artifact. A further evidence in favor of this 
was received by an additional simulation performed for 
H = 4,K = 10,7^ = 10000, which gives very slightly 
shorter core length as H = 4, k ^ 10,7] = 1000. 

We note that, again, our calculation gives slightly 
larger values than those reported for the Faddeev-Skyrme 
Hopfion in [2l[. However, this falls well within the esti- 
mated accuracy of the program and the results show a 
clear approach towards a limiting value, which is more 
important than the exact value of that limit. 

Because every studied soliton follows the same pattern 
of decreasing core length with decreasing k, and increas- 
ing core length with decreasing rj, we conjecture that this 
is a general feature of this model: for core length C one 
has ViJ G Z : lim,,^o C{H) = 00 and lim^^o C{H) = 0. 
If this is true, it raises a tantalizing possibility: start 
with a Faddeev-Skyrme model Hopfion and then begin 
decreasing k and rj in such a way that the collapsing and 
exploding effect of their reduction balance. It is unclear 




(a) ft = 10 (b) r? = 10 

Figure 4: Plots of the normalized core lengths (C/ir3/4) 
of the minimizers for each H = 1 (blue circle), H = 2 
(green triangle), H = 3 (red square), H = A (cyan 
diamond) and H = 7 (magenta pentagon) . We also 
include a single data point for H = 5 (yellow cross) and 
H = 6 (black plus). 



whether such a procedure is possible, but if it is, will it 
give a stable Hopfion at the limit k = 77 = 0? That is, a 
Hopf soliton in pure Ginzburg-Landau model. If such a 
soliton can be constructed, it seems to require an exact 
balance between the growing and shrinking effects of de- 
creasing 77 and At, and as such, is probably not possible 
numerically. 



D. Region supporting stable solitons 

The stable/unstable border is very difficult to ana- 
lyze because the existence of nontrivial local solutions 
depends on the topological invariant, k and i], and their 
detection depends on the initial state. When the start- 
ing point of a simulation is an old solution instead of the 
configuration (|6]) and (O, the stability of the simulation 
also depends on how much k (or r/) is changed from the 
values used to produce the starting point. 

This is best illustrated with charge Q = 4 solutions. 
The minimizers obtained from an initial configuration 
^ using parameters p ~ 1 and g = 4 have a completely 
different shape than those obtained from p = q = 2. 
They also have consistently higher final energies, but we 
were able to find stable solutions for smaller values of 
K, 77 with the first choice of parameters. This reflects 
on the fact observed by Sutcliffe [llj that the number 
of local minima for the Faddeev-Skyrme model increases 
with increasing H . It appears the same occurs in the 
present model as well. It is interesting to note, that the 
most stable soliton is not always the one with the lowest 
energy. 

In order to better understand the relationship between 
the two models, it would be interesting to see, for each 
value of H , which one of the known Hopfions of the 
Faddeev-Skyrme model is most stable in the present 
model. Our single datapoint in this would imply that 
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higher value of E/C would sometimes provide a more sta- 
ble soliton: the configuration p,q = 1,4 has higher E/C 
and can be followed to a lower value of k than p,q = 2,2. 
This would give positive support to the conjecture by 
Babaev that solitons with longer cores would be more 
stable [13. 



V. CONCLUSIONS 

We have studied the modification of Ginzburg-Landau 
model proposed by Ward We find that the stable 

solitons exist for all values of the Hopf invariant H up 
to at least H = 7, but, just like in the situations studied 
earlier [iJ: [3 , this is only possible when the values of 
the parameters k and rj are large enough, and for smaller 
values, the solitons become unstable against Derrick-type 
scaling. 

The results suggest a conjecture that in this model, 
the solitons collapse as k — > 0, but expand without limit 
as 77 — >■ 0. It remains an open question, whether these 
effects could be used to precisely balanced each other and 
provide a way of constructing a stable knot soliton in the 
pure Ginzburg-Landau model by starting from a known 



knot soliton at the Faddeev-Skyrme limit and reducing 
the two parameters until k = 0. 

We gain further insight into the relationship between 
the Ginzburg-Landau and Faddeev-Skyrme models, by 
noting that the order of the values of the normalized 
energies at different values of the Hopf invariant is the 
same as in the Faddeev-Skyrme model [B-0, HH and ex- 
tended Faddeev-Skyrme model [23[. The fact that in the 
modified Ginzburg-Landau model the solitons are local 
minima, not global as in the Faddeev-Skyrme model, is 
interesting from condensed matter physics point of view, 
where local minima are often of great importance. 
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